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Abstract 

We study numerically and theoretically the d-dimensional Hamil- 
tonian motion of fast particles through a field of scatterers, modeled 
by bounded, localized, (time-dependent) potentials, that we refer to 
as (in)elastic non-dissipative Lorentz gases. We illustrate the wide ap- 
plicability of a random walk picture previously developed for a field 
of scatterers with random spatial and/or time-dependence by apply- 
ing it to four other models. First, for a periodic array of spherical 
scatterers in c? > 2, with a smooth (quasi)periodic time-dependence, 
we show Fermi acceleration: the ensemble averaged kinetic energy 
(IIp(OIP) grows as t'^^^. Nevertheless, the mean squared displacement 
(II'Z(^)IP) ^ behaves ballistically. These are the same growth expo- 
nents as for random time-dependent scatterers. Second, we show that 
in the soft elastic and periodic Lorentz gas, where the particles' energy 
is conserved, the motion is diffusive, as in the standard hard Lorentz 
gas, but with a diffusion constant that grows as ||po|| '% rather than only 
as IIpoII- Third, we note the above models can also be viewed as pulsed 
rotors: the latter are therefore unstable in dimension d > 2. Fourth, 
we consider kicked rotors, and prove them, for sufficiently strong kicks, 
to be unstable in all dimensions with ~ t and ~ t^. 

Finally, we analyze the singular case d = 1, where remains 
bounded in time for time-dependent non-random potentials whereas it 
grows at the same rate as above in the random case. 

'Benedicte.AguerQmath.univ-lillel.f 
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1 Introduction 



In this paper, we shall be interested in the classical dynamics generated by 
Hamiltonians of the form 
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H{q,p,t) = ^ + XV{q,t), (1.1) 

where the potential V is periodic in its (dimensionless) spatial variable q 
and periodic or quasi-periodic in its (dimensionless) time variable t, with 
a smooth time-dependence. More specifically, we consider potentials of the 
form 

y{.Q,'t)= W{q - Xm,ujt + (l)o), Xm = '^miei, (1-2) 

where the ei,i = 1 . . . d form a basis of M'^, with ||ei|| = 1 and where W{q, cj)) 
is a function on R'^ x which is spherically symmetric in g, and of compact 
support contained in the ball of radius 1/2 centered at (7 = 0. It is periodic 
in (j) so that, with oj G M"* a frequency vector, the potential is periodic or 
quasi-periodic in time. The model can be thought of as a soft inelastic and 
periodic Lorentz gas, in which the particle with position q and momentum 
p scatters off periodically placed scatterers, modeled by the identical and 
localized time-dependent potentials W . We consider an ensemble of particles 
with fixed given initial energy ||po|P >> ^ starting off in a random initial 
direction from a position close to the origin. We will give numerical and 
analytical evidence for the fact that, for d> 2, the averaged kinetic energy 
(||p(t)P) and mean squared displacement (||Q'(t)|P) of such an ensemble 
behave like 

{\\p{t)f)-t'/^ {\mf)-t'. (1.3) 

The motion is in this sense ballistic, and this in spite of the fact that the 
averaged kinetic energy of the particle increases with time, so that the par- 
ticle undergoes Fermi acceleration. This is, as we will see, a consequence of 
the fact that the particle turns while traveling. We wish to stress that the 
power laws in (|1.3p are identical to the ones valid for random time-dependent 
potentials, studied in [ABLPIO] . In that paper, the motion of a particle in 
a random potential of the form 

V{q,t) = ^XmW{q-Xm,ujt + (j)^) (1.4) 

m 

was considered. Contrary to the situation described in (jl.2p . the scattering 
centers in (|1.4|) . as well as the initial phases (pm and the coupling con- 
stants Am are i.i.d. random variables. It was shown in [ABLPIO] that the 
statistical properties of the motion in such a potential are well described by 
a random walk; this allows for the prediction of the asymptotic behaviour 
of the particles' mean position and momentum. In the numerics performed 
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in [ABLPIO] it was already briefly noted that even if the coupling constants 
and scattering centers are non-random, so that the only randomness remain- 
ing is in the phases (pm, this asymptotic behaviour is unaltered. We go a 
step further in the present paper: if the spatial dimension of the system is 
greater than one, and if the random phases (pm are all chosen equal, leading 
to a purely deterministic potential, periodic in space and (quasi-)periodic 
in time (as in p.2|) ). the asymptotics of the particles' motion is not altered. 
This is, as we will see, a result of the fact that the geometry of the periodic 
lattice together with the instabilities in the dynamics of the individual scat- 
tering events undergone by the particle suffice to effectively randomize the 
motion, even if the potential itself is completely deterministic: the motion 
is in this sense self-randomizing. 

We will support the results in ()1.3p with numerical data for two simple 
models of W on the one hand (Figure [3]) , and with a straightforward exten- 
sion of the ideas developed in [ABLPIO] which, as pointed out already, dealt 
with random time-dependent potentials. The central idea of this analysis is 
the claim that the particle motion is efficiently described by a random walk, 
in which one step corresponds to one crossing by the particle of an effective 
correlation length of the potential, which is present, as we shall see, in spite 
of the fact that the potential is not random. We will show the growth ex- 
ponents in ()1.3p are directly related to the typical momentum transfer Ap 
undergone by the particle during one such crossing. In the present situation, 
it obeys Ap ~ , as we will show. 

In order to illustrate the versatility of our approach we mention it also 
applies to other types of inelastic Lorentz gases that have been considered in 
the literature. In |LRA99l ILRAOOj . the scatterers are hard obstacles, as in 
the usual elastic Lorentz gas, but with periodically or stochastically moving 
boundaries. It is found that the particles' kinetic energy grows linearly in 
time for the stochastic case. This can also be inferred from our analysis, 
since, as the authors show, in these systems, the typical momentum transfer 
in one collision is of order \\p\\^ ~ 1, and does therefore not decrease with 
momentum. As such, these systems are similar to kicked rotors, discussed 
below, for which the energy growth is indeed linear, as we shall establish. For 
periodically moving boundaries, the authors report a slightly faster energy 
growth in the systems they study, presumably due to the persistence of 
correlations; no quantitative analysis is made of the numerically determined 
exponent, however. Let us finally mention that no analysis of the behaviour 
of the particle position is performed in those works, either. We conjecture 
it to be of the same type as in the kicked rotor, namely (||g(t)|P) ~ (see 
below) . 

Our random walk analysis furthermore applies to the soft elastic Lorentz 
gas {d > 2), in which the potentials W in (|1.2|) are independent of time. 
Energy is then conserved and hence the kinetic energy of the particle is 
bounded, and constant in time between collisions. We show that the particle 



3 



then diffuses in space, so that 



{\m\\ 
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~ t. 



(1.5) 



In addition, we show the diffusion constant to grow very fast, as ||po|| ^^'^ 
we illustrate these claims with numerical data (Figure [S]). Note that in 
the usual Lorentz gas, with hard scatterers, the motion has been rigorously 
proved to be diffusive as well |BS81j . but the diffusion constant then scales 
as IIpoII) due to the trivial scaling of the motion with the speed. We are not 
aware of rigorous results on the motion through soft elastic scatterers. 

There is a second interpretation of the results described above, to which 
we now turn. There has been continued interest in the study of time-periodic 
or time-quasiperiodic perturbations of completely integrable Hamiltonian 
systems, both in classical and quantum mechanics. The basic problem is 
the one of their stability. Noting that the Hamiltonian in (jl.ip is invariant 
under the lattice translation group Zei + • • • + Zed one concludes it gen- 
erates a Hamiltonian dynamics on T"^ x W^, where T*^ designates the torus 
ipd _ -j^d^ (Zei -[-••• + Zed). With this change of viewpoint the Hamiltonian 
(jl.ip therefore describes a particle moving on a flat torus, perturbed by a 
scatterer modeled by the time (quasi-)periodic potential W. Such systems 
are referred to as pulsed rotors, and are perturbations of the completely 
integrable system consisting of a free particle on a flat torus. They are 
said to be stable if their kinetic energy is uniformly bounded in time, i.e. 
supt \\p{t)\\ < +00. 

Viewing the Hamiltonian (jl.ip as describing a pulsed rotor, the first 
relation of ()1.3p shows that pulsed rotors are unstable when d > 2. As 
will become clearer below, the potential V constitutes an increasingly small 
perturbation of the free particle, as the particle momentum gets larger. 
The unbounded growth of the energy observed in ()1.3p therefore suggests 
that either all invariant tori of the unperturbed system are destroyed by the 
perturbation, or that the momentum variable diffuses between the remaining 
tori, slowly increasing with time at a rate given in (jl.3p . We have not 
explored any further the complicated mechanisms, possibly related to Arnold 
diffusion, by which this may happen; our results are obtained through a 
much simpler analysis of the momentum transfer in a single scattering event, 
which is sufficient to predict the rate of growth of the diffusion mechanism. 

The preceding analysis does not apply when d = 1, where the motion in 
a random time-dependent potential as in ()1.4p is completely different from 
the one in a potential that is (quasi-)periodic in time as in (|1.2p . Indeed, 
we find for the models we study numerically that, for very long times, (see 
Figure ED 



This is in contrast with the behaviour in a random potential of the type 




(1.6) 
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(fLi|) . where one has (see [ABLPlOj l for d = 1: 



(/(t)>~t§, {q\t))r^tf. (1.7) 

We will explain this difference in Section O Note that (|1.6|) indicates the sys- 
tem to be stable. We will provide two theorems (Theorems 13. II and 13. 2p sup- 
porting this last claim for arbitrary long times using basic time-dependent 
Hamiltonian perturbation theory. More precisely, we show with a suitable 
rescaling of the variables that the system is equivalent to a small, smooth and 
adiabatic perturbation of a simple one-dimensional completely integrable 
system: a particle moving freely on a circle. We then use Lie-Deprit pertur- 
bation theory to show momentum is an adiabatic invariant of the system. 
The result actually holds for arbitrary time-dependent potentials, provided 
they are bounded and sufficiently smooth in time: no assumption needs 
to be made on the (quasi-)periodicity of the time-dependence. Since the 
system is one-dimensional, we have no resonances to contend with. These 
results are not surprising: other one-dimensional non-autonomous systems 
have been studied extensively in the literature over the years, such as the 
Fermi-Ulam model |Ula61j in which a particle moves between a fixed and 
a periodically moving wall, and variations thereof. It is generally expected 
and it has been shown in some cases that sufficiently smooth, periodically 
driven one-dimensional systems are stable: the basic idea is that some KAM 
tori survive the perturbation and create barriers blocking indefinite energy 
growth [LL721 ILL921 ILM88] . Presumably, if one assumes the potential to 
be analytic, similar results can be obtained for the model we study here. 

It is finally instructive to notice the difference in behaviour of the pulsed 
rotors considered above with kicked rotors, in which the time dependence of 
the potentials in (jl.ip is very singular, of the form 

V{q,t) = XY, Ht-nHq), (1.8) 

n 

with V a sufficiently smooth function of its argument. For such systems, it 
is easy to write down the Floquet transformation which gives the evolution 
of the system over one period of the potential: 

^{QiP) = -iP')-, where p'=p—XVv{q), q = q + p ■ (1-9) 

In that case, one finds for A sufficiently large, and independently of the 
dimension d, that (see Figure [8]) 

{\\p{t)f)^t, {\\q{t)f)^t\ (1.10) 

This is again easily understood in terms of the random walk picture we 
develop in Section [21 showing the flexibility of the approach developed in 
[ABLPIO] and here. The main difference with the case of pulsed systems 
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resides in the observation that, whereas for kicked systems the momentum 
change undergone by a particle in one period of the potential is of order 
\\p\\^ ~ 1, independently of the size of the initial momentum of the particle 
and of the dimension d (as is clear from (jl.9p ). this momentum change is of 
order for pulsed systems. Indeed, in one period, the particle traverses 

the unit cell ||p|| times, undergoing each time a momentum change of order 
only, as pointed out above. This fully explains the slower energy 
growth observed in pulsed systems, as well as the slower growth of the mean 
squared displacement. 

We end with some comments on the situation in quantum mechanics and 
the relation of our work with the question of quantum suppression of classical 
chaos. It was observed numerically in [CCIF79] that the one-dimensional 
quantum kicked rotor is stable, even for large coupling (A >> 1), in contrast 
to its classical counterpart, which, as pointed out above, is stable only at 
small coupling. This is the phenomenon that has been dubbed "quantum 
suppression of classical chaos" ; we refer to |Izr93] [Fis93j for reviews on the 
quantum kicked rotor in one dimension. Although stability of the kicked 
quantum rotor for weak coupling was proved in |Bou02j , the problem remains 
open for strong coupling, therefore leaving a rigorous proof of the occurrence 
of quantum suppression of classical chaos in this system open. 

For quantum pulsed rotors, on the other hand, the situation is somewhat 



better understood. First, in one dimension, it follows from [DSSV08| that 
if the potential is small, time-periodic with suitably Diophantine period, 
and sufficiently smooth in both its variables, the system is indeed stable 
|Duc09j . The proof relies on a quantum version of KAM theory and exploits 
the widening energy gap between successive eigenvalues of the unperturbed 
quantum system, typical of the free particle on a one-dimensional torus. It 
is however not clear whether this is an example of "quantum suppression of 
classical chaos", to the extent that, as we argue here, the classical pulsed 
rotor is also stable in one dimension, and this independently of the strength 
of the potential. 

For quantum pulsed rotors of arbitrary dimension, on the other hand, 
various almost stability results have been proved by Bourgain, who showed 
they have a very slow energy growth bounded by t^, for any e > 0, provided 
y is a bounded C°° function of both its variables (see |Bou99aj ). and this, 
whether it is time (quasi-)periodic or not (see also [Del] ) . If the potential 
is analytic, small, and quasi-periodic in time, then in two dimensions this 
growth is even bounded by a suitable power of In t under some additional 
Diophantine conditions on the frequency vector (see j Bou99b| ). The same 
logarithmic bound holds in one dimension, without the smallness hypothesis. 
It furthermore holds for time-periodic potentials in all dimensions, without 
a smallness condition or conditions on the frequency vector. 

These quantum "almost" stability results for pulsed rotors stand in con- 
trast to our results here on their classical counterparts, which we argue to be 
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Figure 1: A particle at time tn impinging with momentum p„ and impact 
parameter bn on the nth scatterer, centered at the point x^^. 



unstable if d > 2, with a power law growth of their kinetic energy given by 
t^/^. As such, these models will provide examples of "quantum suppression 
of classical chaos" provided a fully rigorous proof can be provided of our 
results on the higher dimensional systems considered in this paper. 

The rest of the paper is organized as follows. In Section[2]we describe the 
random walk approximation we use to analyze the motion of the systems 
studied when d > 2. We also show how to derive the various power laws 
for and from it and present our numerical data. In 

Section [31 we study the one dimensional problem. In Section HI we establish 
a comparison with kicked systems. 

Acknowledgments: SDB wishes to thank his regretted friend and col- 
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2 Inelastic and elastic Lorentz gases in d > 2 

In this section, we first briefly explain how the random walk description of 
the motion through an array of random scatterers, developed in [ABLPIO] . 



7 



can be adapted to the present case. We suppose d > 2. We will then sum- 
marily show how the analysis of this random walk leads to the asymptotics 
()1.3p . referring to [ABLPIO] for more details. We consider a typical tra- 
jectory of a particle moving in the potential V described in (|1.2p . Such a 
trajectory will successively "visit" scattering centers X]\f^ at a sequence of 
instants t„, with incoming momenta Pn and impact parameters 6„. Those 
quantities are related by (see Figure [T|) : 

Qn = XN„ - 7:(^n + bn, = jr-^ 6„ • 6^ = 0, ||6„|| < 1/2. 

2 IK II 

Note that we view the impact parameter as a vector perpendicular to the 
incoming direction. We choose to = and xq = 0: the particle starts at the 
scatterer at the origin. We also assume the initial velocity of the particle is 
large, meaning ||po|P >> -^^ At each scattering event the momentum of the 
particle underg oes a change 

= pn+i - Pn- The momentum pn+\ with 
which it leaves the n-th scatterer is also the one with which it impinges on 
the next one. This momentum change depends on the time of arrival tn of 
the particle at the scattering center, via the phase 0^ = <;^o + '^tn; as well as 
on its impact parameter hn- One has 

Ap„ = R{pn, hn, 4'n), (2.1) 

where 

r+oo 

R{p,h,(t)) = -\ drVTy(g(r),a;r + 0) (2.2) 
in which r — )• ^(t) is the solution of 

q(t) = -XVW{q{t),cot + (/.), q{0) = b-\^, p{0) = p. 

2 \\p\\ 

After leaving the support of the nth scatterer, the particle travels a distance 
r]n to the next scatterer, which it reaches after a time Atn = r/n/||Pn+i||- 
The distance hence the time tn+i and the phase (pn+i depend on the 
geometry of the periodic lattice as well as on the dynamics of the scattering 
event via the precise point q^ where the particle leaves the n-th scatterer, 
as well as through the outgoing direction Cn+i- As a result, it is legitimate 
to expect that, after a large number of scattering events, the phases (pn and 
impact parameters 6„ are randomized, that is, they are uniformly distributed 
in their natural range, and that they display short temporal correlations 
only. This hypothesis was tested numerically, and the results are shown 
in Figure [2l In all the numerical results shown in this paper for higher 
dimensional systems, we employed a two-dimensional hexagonal lattice with, 
for N = {Ni,N2) G Z^, xn = N1U + N2V, where u = (1,0), v = ^(l,\/3). 
The potential W was taken in the form of a time-dependent, flat circular 
potential, 

Wiq,t)=fit)x(^^y qeR', (2.3) 
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Figure 2: Distributions and correlations of 6„ and (pn-, with n = 50000, for 
the model described by equation ()2.3p where f{t) = cos(t), calculated with 
10^ trajectories. 



where x{x) = lifO<x<l and x{x) = otherwise. The parameter q^, 
satisfies < < 1/2 to ensure the system has finite horizon and A = 1/6. 
Two different cases were explored for the function /, namely, 

flit) = cos(t) and f2{t) = cos(t) + cos 

In the numerical calculations, each particle has the same initial speed and 
is initially placed at random at a point on the boundary of the scatterer 
at the origin, with an initial velocity drawn with equal probability from all 
possible outward directions. 

The data in Figure [2] strongly indicate that the distributions of phases 
and impact parameters are indeed randomized after a relatively short time, 
and that they decorrelate, despite the fact that the potential is periodic in 
space and time. Under the described circumstances, the statistical behaviour 
of the ensemble of particles is expected to be the same as that of an ensemble 
executing the random walk 

Pn+l =Pn + R{Pn, K, 4>n), Wl = + T: 9n+l = Qn + f?*en+l, (2.4) 

\\Pn+l\\ 

where 77=,, is the average distance traveled between scattering events, g„ the 
position of the particle at the nth such event and in which bn and (j)n are 
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Figure 3: Numerically determined values of , and 

(Ik(OlP) ™ a two-dimensional hexagonal lattice for the potential described 
by equation ()2.3p where f{t) = cos(t). In each plot, different symbols corre- 
spond to different initial conditions, as indicated. 



treated as independent random variables, each distributed in its natural 
range. We will write (•) for the corresponding averages. 

The analysis of this random walk was performed in detail in [ABLPIO] 
and in [AgulO] , We reproduce here only the essence of the arguments al- 
lowing to predict the asymptotic behaviour of the particles' averaged energy 
and position, referring to the above works for more mathematical details. 
We will in particular only carry the leading order terms in all expansions. 

We first describe the growth of the kinetic energy of the particle (see 
Figure [3l top panels) . Leading order perturbation theory applied to ()2.ip 
and ([22]) yields 

IK+ill ~ IKII + f^, where/3W = A / d^ (uj-V^)W{b + 

\\Pn\\ J — oo 

(2.5) 

Since, if the (pn are uniformly distributed on the torus, {f^n'^) = 0, one can 
think of (j2.5p as a random walk with variable step size and it is then easy 
to determine the asymptotic behaviour of ||pn||- Indeed, since 

A||p„f ~3||Pnf A||p„|| ~3/3«, (2.6) 

executes a standard random walk with ||p„||-independent step, so that 

(l|Pnf)~n^, and (||p„|n ~ nf (s > -3), (2.7) 
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since the correlations between the Pn '^ are supposed to decrease fast. Com- 
bining this result with (|2.4p . one finds 

(t„)~ni, so that (||p(t)|n~tt, (2.8) 



which proves the first equation in (jl.3p . This slow growth of the energy 
(s = 2) is indeed observed in the numerical data presented in Figure [3] (top 
panels). It is identical to the one found for random time-dependent poten- 
tials in |ABLP10] confirming that the motion in a time-periodic potential of 
the type considered is self-randomizing due to the instabilities inherent in 
the scattering mechanism. 

To understand the asymptotic behaviour of the particle's position q{t) 
we have to analyze how fast it turns. For that purpose, notice first that 
()2.2p and ()2.4p allow one to write, to lowest order in perturbation theory: 

r+oo 

Pn+1 '^Pn+ ip-n") where a^^) = "-^ / (^n + f^^ri, 4>n) ■ (2.9) 

\\Pn\\ J — oo 

Note that a^n^ • e„ = and (ai^^) = 0. It follows from ([23]) and i^Mj that 

en+1 ^ e„ + rr^T. (2.10) 

\\Pn\r 

One can therefore think of the e„ as executing a random walk on the unit 
sphere in R'^. For values of m small enough that ~ \\Pn\\^ that is for 

m « n (see (|2.6p - (|2.7p ). it is legitimate to approximate this walk by 



m— 1 



k=0 

Thinking of the a^+fc independent steps (again, for a more careful treat- 
ment, see again [ABLPIO] ). it follows that (||en+m — ~ m^/n^. We 

2 

conclude that (11 

(^n+m ^n||) ~ 1 fd" times of order ?n, ~ ?7.3 << re. This 
means in particular that, if after n collisions, we consider only those parti- 
cles moving in a given direction e, so that en = e, it will take of the order of 
7^2/3 jj^Qj.g collisions for the distribution of particle directions to spread over 
a macroscopic angle about e. This spreading is qualitatively illustrated in 
Figure m where we see that after m ~ 1500 collisions, the set of directions 
{sk} keli-mj covers the unit circle. 

Writing for the number of collisions needed for the particle to turn 
k times over a macroscopic angle, we find 
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(d) m = 1500 

Figure 4: Coordinates of the vectors e^, k G [1; m|, with m as indicated, for a 
single trajectory in the model described by equation (j2.3p with /(t) = cos(t). 
The initial velocity is taken to be po = 2eo, and eo is represented by the large 
dot. 



so that Mk ~ A;^. It is then reasonable to think of the motion as being 
effectively rectilinear at constant speed between and M^+i, so that one 
obtains an effective random walk for the particle position given by 

QMk+i = QMk + f?*(Mfc+i - Mk)eMk - QMk + il^k^euk- 

Treating the e Mfe as independent random variables distributed uniformly on 
the sphere, one finds from this that 



Extrapolating to all n, and using (j2.8p yields the second equation in (|1.3p 
(see Figure O bottom panels). 

We point out that the previous analysis about the evolution of g„ applies 
equally well to a time-independent potential W . In that case, the kinetic 

1 1 2 

energy ||pn|| of the particle is of course a constant of the motion and one 
finds, instead of (l2lil 
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Figure 5: Numerically determined values of (left panel) and of 

the diffusion constant (right panel) in the case of a soft Lorentz gas with a 
potential as in ()2.3p . with /(t) = 1 and A = 0.49^/2, for 12 initial velocities 
from 0.5 to 3.8. 



so that 

1 '""^ A 

en+m = Cn + — TT^ a^+fc, and (||en+m - e„||) ~ - — [T^\/rn, 
llpoir ^ llfoir 

provided ||po|| is large, that is ||po|P >> ^- The directions e„ now diffuse 
on the unit sphere at a fixed rate, that decreases with the particle's speed 
as IIpoII"^: indeed, the faster the particle, the less it turns in one scattering 
event. It follows that Mfc+i = + A'^Hpoll"^, so that ~ A~^ || po T ^ 
from which one concludes with a reasoning as above that 




In other words, our random walk analysis shows that in a soft elastic Lorentz 
gas, the particle diffuses, and that the corresponding diffusion constant 
grows like ||po||^- This is illustrated in Figure El One notices in the left 
panel that at high speeds the diffusive regime sets in at later times, as ex- 
pected: this is apparent from the bending of the topmost graphs. For short 
times, the motion is ballistic, of course. It is furthermore interesting to note 
the difference with what happens in the hard Lorentz gas, where diffusion of 
the particle has been rigorously proven; one has ||g(t)||^ ~ IIpoP- The slow 
growth of the diffusion constant as ||po|| can in that case be explained by 
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the observation that, in the presence of hard obstacles, the trajectories do 
not depend on the speed of the particle, so that the velocity dependence of 
the diffusion constant follows from a simple scaling argument. 

This ends our analysis of the class of soft non-dissipative inelastic and 
elastic Lorentz gases in dimension d > 2, described by ()1.2p . As pointed 
out in the introduction, the Hamiltonians ()l.ip - ()1.2p can also be viewed as 
describing pulsed rotors: our results here show they are unstable in higher 
dimension since their energy grows on average. A comparison with kicked 
rotors will be made in Section HI But we first turn in the next section to 
the study of the one-dimensional situation, which is very different. 

3 Inelastic Lorentz gases and pulsed rotors in one 
dimension 

Unlike the standard (hard) Lorentz gas, which has no interesting analog in 
one dimension, the soft Lorentz gas described by ()l.ip - ()1.2p yields a mean- 
ingful and non-trivial dynamics in one dimension as well. It is to its study 
we turn our attention in this section. 

We have just shown that, in dimensions strictly higher than one, the 
asymptotic behaviour of a high speed particle in a potential periodic in 
space and (quasi-)periodic in time (as in ()1.2p ) is identical to its behaviour 
in a field of random scatterers (as in (|1.4p ). We will now see that in one 
dimension the situation is totally different. It was shown in [ABLPlOj that 
in the random case, the particle momentum grows as , just as in higher 
dimensions. We will show here that, on the contrary, the momentum of a 
particle in a space and time periodic potential remains bounded for very 
long times. The precise results are given in Theorems 13. II and 13. 21 below. In 
fact, our conclusions hold for spatially periodic potentials with an arbitrary 
smooth and bounded time-dependence: time-periodicity is not needed. The 
results follow from time-dependent Hamiltonian perturbation theory, with 
the inverse of the particle's initial speed as the small parameter. 

To understand the notion of high speed in the present context it is helpful 
to remember how the dimensionless Hamiltonian H, position q, momentum 
p and time t in (jl.ip arise. The problems we are considering have a natural 
length scale i, which is the size of the unit cell of the lattice, and a natural 
time scale T, which is the typical scale on which the potential fluctuates. 
This suggests expressing the position q of the particle as a multiple of i, its 
momentum as a multiple of i/T (we set the mass equal to 1) and to measure 
time in units of T. So, writing q = iq, p = pi/T and r = tT, the equations 
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of motion for {q,p), equivalent to those obtained from (|l.ip are 




(3.1) 



dp, , Xi 




( 



I 'T 



) 



which are Hamilton's equations for 



t T 



q T 



) 



with = A£^/T^ having the dimension of an energy. In what follows, 
the potential y is a dimensionless smooth function in C°°(M^), periodic of 
period 1 in its first (spatial) variable, and bounded with uniformly bounded 
derivatives in both its variables. It need not be of the specific form given in 
p.2p and in particular, we do not assume the potential is (quasi-)periodic 
in time. As before, the spatial periodicity of the potential implies we can 
equally well think of the equations as describing a particle moving on a circle 
of circumference £ in a bounded time-dependent potential. 

This Hamiltonian has one dimensionless parameter, which is A = /xT^/-£^, 
where /i is the typical strength of the potential and ^^/T^ the kinetic energy 
of a particle crossing the unit cell in a time T. We will study the solutions of 
the equations of motion for large initial momentum p(0), by which we mean 



for some fixed 2 < a < A < +00 and for a large value of p.^, to be specified 
below; is therefore a second dimensionless parameter, measuring the rel- 
ative size of the initial momentum p(0) and i/T. Note that e = l/p^t is the 
time a free particle with initial momentum p^ needs to traverse the unit cell 
of length 1, so that e << 1 means this time is short compared to the time 
interval 1, characteristic of the variations of the potential. The particle is 
in this sense fast. We will in addition suppose Ae^ << 1, which means that 
the kinetic energy of the particle is much larger than the typical size of the 
potential: so the particle is energetic. 

What will happen is intuitively clear: since the particle is energetic and 
fast, the potential provides a small and slow perturbation of the free motion. 
In fact, in one traversal of the unit cell, the particle will pick up a momentum 
change of order p~^ at most (see (j2.9p ). Naively, one may expect that 
after a time of order 1, the particle has turned around the circle p^, times, 
resulting in a momentum change of order 1 . If the momentum changes were 
systematic, this could lead to an acceleration p{t) ~ t, much faster than 
in a random time dependent potential where it was proven in [ABLPIO] 
that p{t) ~ tV5. As we shall prove, the situation is much better than that: 



a'P*7^ < p(0) < Ap^— 4^ ap^ < p{0) < Ap^ < +00 
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Figure 6: Numerically determined values of and q{t)'^ in one dimension 
for the potential V{q,t) = cos tX[o,2/3] Iii each plot, different symbols 
correspond to different initial conditions, as indicated. 

due to cancellations, the momentum is bounded for arbitrarily long times 
(Theorem 13. ip . The phenomenon is illustrated in Figure [6] for the potential 
V{q,t) = costX[o,2/3](9)- Note that this potential is not even smooth in q, 
which doe s not seem to affect the stability of the momentum at all, on the 
time-scale explored. The difference with the situation in higher dimensions is 
intuitively clear: in one dimension, there are no impact parameters that can 
be randomized by the successive scattering events and in addition, since the 
distances between successive scattering events are all rigorously identical, 
the phases of the potential do not randomize either. 

We now give precise statements of the above claims, supposing 

Vt G R, / V{q,t) dq = 0. (3.2) 
Jo 

This condition constitutes no restriction since we can always substract its 
spatial average from the potential without changing the equations of motion. 

Theorem 3.1. Let 2 < a < A < +oo. For all a^, > and all K eW, there 
exists a constant Ck > such that, for all large enough and for all initial 
conditions ap^ < p(0) < Ap^,,q{0) G [0, 1[, one has, for <t < <J^p^~^, 

\p{t)-pm<CK ^ 



p{0)- 



The error estimate is optimal in this result, as Theorem 13. 21 below shows. 
The variations in p(t) do really have an amplitude of ^(0)~^: 

Theorem 3.2. For all a^, > and all < t < ct*, p{t) = Pap{t) + O (p^^) , 
where pap{t) = p{0) - A (F (g(0) + p{0)t, t) - V (g(0), 0)) /p(0). 
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(a) Evolution of p{t). 



(b) p{t) VS Pap{t). 



Figure 7: On the left, the evolution of p{t) for ten different initial velocities, 
from 0.2 to 2, for the potential V{q,t) = costX[o,2/3]{Q)- On the right, 
numerically determined values of p{t) compared to values of Pap{t) for the 
same model and with initial velocity p{0) = 2. 



This result is illustrated in Figure [71 Its proof provides in addition an ap- 
proximation of p{t) to order for arbitrary long times (See Corollarv l3.3p . 
We will prove these results using standard techniques of Hamiltonian per- 
turbation theory as explained in [CarSH ILL831 ILM88j , only sketching the 
essence of the arguments and stressing the points particular to the present 
model. Details are given in lAgulO]. We note that Theorem 13.11 is close 
to a statement that can be found on p200 of |LM88]. We first rewrite the 
equations of motion in a form suitable for Hamiltonian perturbation theory. 
For that purpose, we define a new time-scale and new variables 

a = - = p^t q{a) = q{ea), p{a) = ^^if^ = ep{ea). (3.3) 
e P* 

Here a counts the number of revolutions executed by a particle of momentum 
in a time t. Then, if {q{t),p{t)) solves Hamilton's equations corresponding 
to (jl.ip . they satisfy 



^(^)=P», 

^^a) = -e'XdiV{q{a),ea), 



(3.4) 



with initial conditions q{0) £ [0, 1[ et 2 < a < p(0) < A < +oo. These are 
Hamilton's equations corresponding to H defined by 

H{q,p, e, r) = ^p^ + e^XV{q, r), (3.5) 
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where r = ea. We will write {q{a),p{a)) := {q{0) , p{0)) for the Hamilto- 
nian flow generated by H. 

Proof of Theorem \3.1[ One considers a C°° function x{QiPj ''") periodic in q, 
defined for p > 0, fi,T G M, which is bounded together with all its derivatives 
on p > 5, for all 6 > 0. We will show a judicious choice of x allows to use 
its flow to define new coordinates {Q, P) in which the equations ()3.4p take 
a simple form. More precisely, we define the local diffeomorphisms 

$J^^^)(g,P) = (g-(/i,T),p(^,T)), 

through the Hamiltonian flow in ^ of 

"^^[f ' = {dpx) {q{l^,T),p{fl,T),fl,T), 

— — = - (dgx) iq{tJ-,T),p{n,T),n,T), 

with initial conditions (g(0, r),p(0, r)) = {Q,P). If /i is small enough, then 
for all T and for k = 1,2, 

^ ■^>^ [a/k;oo[-^Rx [a/2k;oo[. (3.7) 

At the end of the proof we will take = e, which means p^, = 1/e has to be 
large enough. For all /i < e and all a small enough so that {[a, +oo[) C 
[a/2,+oo[, we define for (Q,T) e M x [a/2,+oo[: 

We also define 



(Q(0),P(0)) = ^ci>x^oJ (g-(0),p-(0))and(g(a),P(a)) = V'(^,,)(Q(0),P(0)). 

Since V'(^,ct) is symplectic, there exists a Hamiltonian H{Q, P, fi, r) of which 
V'(/i,(T) is the flow with respect to a, with ^ fixed: 

^/ ° = {/; h] o V/ : M X m+ ^ m. 

We will write V'(^,(t) = • The expression for is known ( |Car8H ILM88) ): 

^ = ^°^f,,.)- r^'^Xo<l>5,,,)d^'. (3.8) 
We obtain the result by constructing x as a polynomial of degree K in n, 

K 

x{Q,P,fJ','r) = ^/i"xn+i(g,P,r), 

n=0 
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for each K £N* , in such a way that for fJ- = £, 

^P{a) = -dQH{Qia),Pia),e,T) = 0(e^+i); sup |P(a) -p{a)\ = 0{e^). 

(3.9) 

To obtain this, the smooth functions Xn wih be chosen so that, when /U = e, 
H is independent of Q up to terms of order e^^^. Clearly, for all / G 
C°^(M X M+,M) and for all r > 0, g, Q G P > 0, the functions 

^ ^ / o $ J^^^^ (Q, P) and fi^ fo J.^^)) ^ (g, p) 
are C°°. Indeed, one easily sees that, with := {•; x} + d^, 



A similar argument works for /i i— )• / o |^<I>^^^jj {q,p)- Hence, such func- 
tions can be expanded to any order K, with smooth coefficients in {q,p,T). 

We now take = e in I^Bj- With H as in ([33]), ^ is in e. We 
introduce the notation 



K 



H = Ho + e^H2, if = ^e"i?„ + 0(e^+i) = /7^ + 0(e^+i), (3.10) 



n=0 

K 

n=0 

1 



n=0 



Taking the derivative of (|3.8p with respect to e and composing with ^) 

on the right, one finds daX = f^e-f^ + {^^x} ~ deH o ^<I>^^^^^ . Since 
da = £dr, one then easily identifies the terms in e" for n € |0;i^ — 1], on 
both sides of this equation. For n = 0, one finds dgXi = —Hi/P, since 
Hi = 0, Hq = p^ /2 (see ()3.5p ). Since we wish to choose Hi to be indepen- 
dent of q, averaging this equation over q implies Hi = 0, and hence dqXi = 0; 
one chooses the particular solution xi = 0- For n G [1; X — 1], one has 

n-l 



drXn = {n+l)Hn+l + {HQ-Xn+l] +^{Hn~k;Xk+l] 

-(n + l)ToHn+i - (^) - ^ + ^)TkH., 
k=i ^ ^ 



n—k+l ; 
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where = when k ^ {0; 2}. Since Tq = Id one finds, 

n-l 

(n + 1) (Hn+l - Hn+lj - drXn + ^ {Hn-k;Xk+l} 



dqX 



n+1 



1 

P 



k=0 



E 

k=l 



{n-k + l)TkHn-k+i 



(3.11) 



For k G [l;ra], the coefficients Tkf depend only on the Xj with j £ [1;^]. 
Indeed, we aheady know Tq = Id and, in addition, the equahty 



;x 



impHes, after identification of the terms of order /i" for n G N, 

n 

(n + i)r„+i/ = -J2 {Tn-kf; Xk+i} ■ 

k=0 

After averaging over q in (|3.11|) . one can therefore recursively construct the 
coefficients Xn in such a way that the coefficients Hn,n G [0;K], don't 
depend on Q. This yields the first part of (|3.9|) . 

To show the second part of ()3.9p . we note that, since xi = 0, (|3.6p . ()3.7p 
and the definition of P{cr) imply 



p{a)-Pia) 



Now, ()3.1ip for n = 1 reads dgX2 
grating this equality with respect to q yields 



r {Ii^gX2{q{^^,r),p{^l,T),T) + 0{fi^)) dfi. (3.12) 
Jo 



2 Ho -Ho 



/p, since xi = 0. Inte- 
-^H2 since H2 has zero 



average and H2 is to be chosen independent of q. As a result, 

dgX2 = bl2 = 

p p 

The hypotheses on V together with ()3.7p therefore imply sup^>o |p(o") — -P(o") | 
0{e^). Similarly, sup^>Q [^(cr) — Q((t)| — ©(e'^). One concludes the proof of 
the theorem by coming back to the definition of p (see (|3.3p ): 

Pit) - piO) = p, ipia) - pm = e-^ ipia) - p(0)) . 

This difference can be decomposed in three terms 

\p{cT)-pm < \p{a)-P{a)\ + \P{a)-P{0)\ + |P(0) -p(0)| . 



Using (|3.9p . one concludes there exist constants ck et Ck such that 

\pit)-pm < e~^CK{e^ + e'^+^a + e^)<CKe (3.13) 
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for all < t = ecr < cr*e--^+2_ 

□ 

Proof of Theorem \3.2[ In view of the definitions of q, p and a (see (|3.3p ). it 
suffices to show that, for all < t = eo" < a*, one has 

p{a) = p{0) - {V im + p{0)a, r) - y(g(0), 0)) + O(e^). (3.14) 

l^et K en, K> 3. From (fXH]) . using P{a) - P(0) = 0{e^+^a), 

p{a)-p{{)) = p{a)-P{a) + P{a)-P{0) + P{0)-p{0) 

-,,^^m/^,r)r)_y(,-(,,0)0)X ^ 



+0(e3) + 0(e^+V). (3.15) 

Here ^(0,r) = Q{a) andp(0,r) = P{a). In addition, since xi = 0, it follows 
from ()3.7p for all a and all < /U < e that 



g((T)-(3(cr) = / apx(9(/^,T),p(^,r),^,r)d/i = O(e^), 

JO 

Q(a)-g(^,r) = 0{e^), p{a) - p{fi,T) = 0{e^). 
Inserting this in (|3.15p . and using p{a) — p(0) = O(e^) (see ()3.13p ). one finds 

= J, (mihii _ mm.) + o(e-) + o(e-«,) 

V pI") Pio) J ^ ' ^ ' 

= (V(Q(''),T) - V(f(0),0)) + 0(e3) + o(.-*-+V). 

Finally, the definition of (equation ()3.10p ) implies 

Q{a) = Q(0)+ r dpH^(P{a'),e,ea')da' + 0{e^+'a) 
Jo 

= Q(0)+ r dpH^(PiO),£,ea')da' + 0{e^+^a) 
Jo 

= g(0)+ r dpH^{p{0),e,ea')da' + 0{e^) + 0{e^+'a). 
Jo 

So 

Pi'T) = PiO) + £ dpH^iPiO),e,ea')da',r'^ - y(g-(0),0)) 

+0(e3) + 0(e^+V). (3.16) 
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To conclude the proof, one fixes K = 3. By definition 

H^ = Ho + eHi + e^H2 + e=^#3, 

and we already know Hq = Hq = /2 and Hi = H2 = 0. Considering p. lip 
with n = 2 yields dqXz = j Sffa — drX2^ ■ Averaging over q implies Hy, = 
0, since X2 can be chosen of zero average. Equation (j3.16p then yields (|3.14p . 

□ 

Equation (j3.16p implies the following corollary. 

Corollary 3.3. Under the hypothesis of Theorem \3.1l for all K N, K > 3, 

and all < t < a^p^~^, 

where 

pfpW =p(o)-^ (v ((?(o) +p,J\hk (^PM^p-\t'^ dt',t) -y(g(o),o)) 
4 Kicked rotors 

We now turn to the kicked system described in ()1.8p - ()1.9p and consider 
{qn,Pn) = ^"(go,Po) G T*^ X so that 



Pn+l =Pn- XVv{qn), 
Qn+1 = Qn+ Pn+l- 



(4.1) 



Note that here g„ = q{n) and pn = p{n) since in (jl.Sp the potential is of 
period 1. So in kicked systems, there is exactly one "scattering event" per 
unit time, independent of the particle's speed, whereas in pulsed systems 
the number of kicks per unit time increases linearly with \\p\\ as the particle 
speeds up. So in kicked systems, the momentum change Akick||p|| is of order 
A in a time interval of order 1, independently of the size of \\p\\. On the other 
hand, in pulsed systems of dimension d > 1, as a result of (|2.5p . supposing 
the are i.i.d., the typical momentum change Apuis||p|| of a fast particle 
with momentum ||»|| in a time interval of order 1 is 




So Apuisll^ll ^ Akick||p||- We will now see that, as a result of this, kicked 
rotors speed up much faster than pulsed ones. From (|4.ip we find 



n-l 



:po-Aj](vy)(g,). 



t=0 
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Figure 8: Numerical determined values of and ^11(7(^)11^^ in one 

dimension (top) and two dimensions (bottom) for the model with kicked 
rotors described by equation ()1.9p . where V{q) = Y[i=i "^^^ ^^'^ ^ — 10- 
The different symbols correspond to different initial conditions, as indicated. 



Since in kicked systems Ap„ is of order A, independently of the size of ||p||, 
the second equation in (|4.1|) suggests that, for large A, these successive large 
momentum changes will lead to a randomization of the position of the system 
on the torus, allowing one to think of the Qn as uniformly distributed on the 
torus, and with short temporal correlations. This implies {{pn ~ Po)'^) ~ n 
or equivalently (||p(i)|P) ~ t. In dimension d = 1, this immediately implies 

^lk(^)ll^^ ~ dimension d > 1, the change of direction e^+i — e^, with 

en ■= Pnl \\Pn\\ , is of Order \\pn\r^ ~ n"^/2, so that (||e„+m - e„||) ~ 1 for 
m ^ n. Writing, as in section [21 Mk for the number of collisions needed by 
the particle to turn k times over a macroscopic angle, we find M^, ~ 2^ and 

(iMk+i = QMk + {Mk+1 - Mk) PMk 
3/2 

This yields (H^MfeP) ~ M| and, interpolating between the M^, we therefore 
finally find (|| (7(^)11^) ~ even in dimension d > 2. The reason for this is 
clear: the particle speeds up considerably faster for kicked systems (where 
< l|p(OlP >~ than for pulsed ones (where < ||p(t)|P >~ t^^^) and hence 
turns much less. These asymptotics are in agreement with the results of 
numerical computations displayed in Figure [HI 
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